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(-H ' Abstract. In this article we create a new algorithm for the perfect simulation 

' of the infinite Potts model at a sufficiently small or at a sufficiently high 

temperature, in particular under the transition phase temperature. We study 
the model for free boundary conditions and we give some consequences for the 
constant boundary conditions. 

Qh . Given a countable graph G = (V, E), a positive number q and parameters p — 

{p e £ [0, 1] : e £ E}, the random cluster measure is defined on the measurable space 
(il, J 7 ), where fi = {0,1} E and T is the cr-algebra generated by finite cylinders. 
This measure was introduced by Fortuin and Kasteleyn as a way to study the Ising 
and Potts models (see [FK72]). 

In our paper we do not require the parameters p e to be all equal to the same 
constant; however we limit our attention only on the models with q > 1. This choice 
has been made both to maintain the article simpler and also because the case q > 1 
has important connections with the statistical mechanics and in particular with the 
models of Ising and Potts. 

The aim of this paper is to construct an algorithm which gives a perfect simu- 
lation of a random cluster measure on a finite region of an infinite graph. Notice 
that even if the perfect simulation is obtained only on a finite region, it takes into 
fT") ■ account the fact that the random field on this region is influenced by the value of 

the filed on the whole infinite graph. On any finite graph it is possible to define, in 
an explicit and computable way, the random cluster measure. However, the random 
cluster measure associated to the finite region in this way, in general, is not the 
restriction of the measure of the infinite graph on this finite region. 
, Now we briefly explain how this simulation is obtained. As will be recalled 

in Sections 2 and 3 the random cluster measure is invariant under a markovian 
dynamics. Introduce a countable number of copies of the graph G and think them 
as placed at level, 0, —1, —2, etc. Choose also an order of the edges of G: e\, e2, 
etc. For a configuration w S £1 of the graph at level — M create new configurations 
at level above — M updating the value of uo £k one at the time, according to the 
conditional probabilities that depends geometrical shape of the configuration. The 
details of this dynamics are given in Section 3. Let us just recall here that the 
law used to update the value of uj ek depends on the existence of a connected path 
of edges e, different from efc, such that ui e = 1, joining the end vertexes of e&. 
The construction of this dynamics can be seen as a particular case of the Glauber 
dynamics. In the study of the random cluster measures similar dynamics where 
already considered, for example by Grimmett in [Gri95] . 

We now construct a coupling of all the dynamics for each possible initial config- 
uration. We color with black, gray or white, in an independent way, the edges of all 
the copies of G. The law used to color the edge e at level —I depends only on the 
parameter p e . The coupling is constructed in such a way that for all possible initial 
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configurations uj, the value of the configuration in e at level —I will be 1 whenever 
the edge e at level —I is black and whenever it is white. 

Now fix a finite region F of G. In Theorem 8 in the case of high temperature 
(which corresponds to p e close to zero), we prove that, for almost all coloring as 
above, there exists a finite region Cp of the union of all the copies of G that is 
"surrounded" by white edges and containing the region F at level 0. Finally, in 
Theorem 5 we prove that, given a coloring, if the region C b F is finite then we can 
determine a bigger finite region H such that the output of the dynamics described 
above at level in the region F does not depend on the choice of oj and on the 
coloring outside the region H. 

Similar results are proved in Theorems 7 and 4 in the case of low temperature 
which corresponds to p e close to one. However in this case the meaning of the 
word "surrounded" is different. To treat this case we have to make some further 
assumption on the geometry of the graph G. For this reason in Section 4 we 
introduce the concept of simplicial graph. These are the graphs that can be obtained 
as the vertexes and edges of a tessellation of an euclidean space. The first example 
we have in mind is the cubic lattice L d = (Z d ,Ed). In Section 4 we prove also the 
results required in the proof of Theorems 7 and 4. 

In Section 7 we summarize these results and we show how to obtain an algorithm 
for the perfect simulation of the random cluster measure. We give the details only 
in the case of low temperature which is the most difficult and interesting. The case 
of high temperature can be treated in a similar way. As a byproduct we also obtain 
the uniqueness of the random cluster measure at low or high temperature, which, 
al least in the case of L d is well known (see [Gri95, Gri04]). 

In the last section we analyze the connection between the perfect simulation of 
the random cluster measure and the Potts model. We construct an algorithm to 
simulate the Potts model with free boundary conditions at low or high temperatures. 

We also refer to some literatures on perfect simulation. In [PW96] was introduced 
the perfect simulation algorithm for Markov chains. If the chain is ergodic with 
this algorithm it is possible to simulate the unique stationary measure associated 
with the chain. This paper has started some new research fields. One area of 
research concerns the Markov fields (see [HS00, DSP08]); a second one concerns 
the processes with infinite memory (see [CFF02, DSP12, Galll]). Recently, these 
two areas of research have been in some sense unified by studying Gibbs measures 
with infinite interaction range (see [GLO10, DSL12]). Our paper is included in the 
latter context. 



1. Some notations on graphs 

In this section we recall some definitions on graphs that will be used in the sequel. 

In this paper a graph will be a collection of two sets, V called the set of vertexes 
and E called the set of edges, and of a map from E to the set of unordered pairs of 
different elements of V. The pair associated to an edge e are called the end vertexes 
of e and the two vertexes arc said to be adjacent. As it is common in the literature 
we will denote a graph by (V,E). A path in G joining the vertexes u and v is a 
sequence ei, . . . , e m in E such that e% and ej+i have a common vertex, u is an end 
vertex of ei and v is an end vertex of e m . The integer m is called the length of 
the path. Two vertexes are said to be in the same connected component if there 
is a path joining them. The graph-distance of two vertexes u and v, denote with 
5g{u, v), is the length of a minimal path of joining them, and it is infinite if the two 
vertexes arc in different connected components. We denote by E>g(v, r) the ball of 
center v and ray r with respect to this distance. 
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2. The random cluster measure 

In this section we define the random cluster measure introduced by Fortuin and 
Kasteleyn as explained in the book of Grimmctt [Gri04] . Since our setting will be 
slightly more general than the one exposed by Grimmett we give the construction 
of the measure. However all the arguments given in his book easily generalize to 
our setting, so we refer to [Gri04], Chapter 4, for the details. 

2.1. Construction as thermodynamic limit. For our constructions we fix a 
graph G = (V, E). We further assume that it is countable of finite degree, meaning 
that V is countable and that every vertex is an end vertex of a finite number of 
edges. 

Set ft = {0,1} B and let F be the cr-algebra generated by finite cylinders. If 
co £ £1 we denote by E{uj) the set of the elements e £ E such that u> e = 1. To 
define a random cluster measure we also fix parameters p = (p e £ [0, 1] : e £ E), 
and q £ (0, oo). For simplicity in this paper we will assume q > 1 which is the more 
significant for the application to statistical mechanics. 

There are two ways of defining random cluster measure on G. The first method 
is as limit on finite subgraph and is called the thermodynamic limit. The sec- 
ond method is by giving the conditional probabilities on finite subgraph and it is 
called the Dobrushin-Lanford-Ruclle or DLR method. We now explain briefly the 
construction as thermodynamic limit. 

Given £ £ Cl and F C E a finite set, let flp = {w £ O : oj e = £ e for all e £ F}. 
We define the measure 4>% p q on O by: 

4 M = (^7[n ee ^e e (l-fe) 1 --]^) if W en|, (1) 

'' 1 otherwise, 

where k(uj, F) is the number of connected components of the graph (V, E(uj)) that 
intersects F and Z^f is just the normalizing constant. 

Following [Gri04] Definition 4.15 we say that a probability <j) on F) is a limit 
random cluster measure if there is a sequence (£„, F n ) such that <fi is the weak limit 
of the measures <ffip and we denote by W p . q the set of these measures. 

If we fix to be constantly equal to 1 (resp. to 0) the limit of the measures (j^p 
exists, and it does not depend on the choice of the sequence F n (see Theorem 4.19 in 
[Gri04]). This limit will be denoted by (f) p q (resp. <pp tQ )- Moreover, for all <j) S W p>g , 



we have 



6p,q <st 4> <st (t>p, q - 



where < st is the standard stochastic order as in [Gri04] Section 2.1. 

Another important property of these measures is the so called finite energy prop- 
erty. Let 

Pe =Pe/(Pe + 9(l ~ Pe)), (2) 

then for all 4> & W p , q and e £ E, we have 

Pe < <t>(L e \T e )(u) < p e , U.S. 

where L e = {uj : uj e = 1}, T e is the cr-algebra generated by the finite cylinders with 
base contained in E \ {e}. 

2.2. DLR construction. In the case of a finite graph G — (V,E) the definition 
given above furnishes a unique measure on (Q, F) which is characterized by the 
conditional probability of u> e — 1 given the values of u> in E \ {e}. In this model 
these probabilities depend on the existence of a path in E(co) \ {e} joining the 
two end vertexes of the edge e. Let K e be the set of configurations ui having this 
property. 
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In the case of infinite graph this property can be formalized as follows: 

^.|T.)M = (* (3) 

\p e if u g K e . 

A probability <f> on (O, J 7 ) is a DLR random cluster measure if it satisfies equation 
(3) for all e E E. We denote by 1Z Piq the set of these measures. 

In the infinite setting the two definitions of random cluster measure are not 
always equivalent. However (see [Gri04], Chapter 4, Section 4) it is known that 
lZp, q is not empty, in particular </>° and (j) pq are elements of 1Z Pyq and for all 



G Tl Pl q one has 



<std><, - 



In particular notice that 

card(>Vp jg ) = card(7£ Pi q) = 1 if and only if <fp p q = 4> p q . 



3. Construction of the dynamics 

In this paper, following the literature, to give a sample of a measure on (O, T) 
we introduce families of auxiliary random variables (u e ) e ^E that are independent 
and uniformly distributed on [0, 1], which in the algorithm we present in Sections 
7 and 8 can be thought as the output of a pseudorandom function on a computer. 

We now define a stochastic dynamic such that the measures in H p , q are invariant. 
A similar dynamics was already considered by Grimmett in [Gri95] . 

Let G = (V, E) be a countable graph of finite degree and we choose an order for 
its edges so that E = {ei,e 2 , . . .}. 

For a negative number N we define An = {(n, k) eZxZ: —1 > n > N and 
k > 1} (rcsp. A = \J N An) and U = [0, 1] A . On U we put the Lebesgue product 
measure so that the coordinates u n .k are i.i.d. random variables having uniform 
distribution on [0,1]. We define also An = An U {(n, 0) : n = N,...,0} and 
A = \J N A N . 

For a fixed N < and for a fixed Xn,o ■ U — > Vl let (X n _k '■ U — > il such that 
(n, k) G An) be a process with values in Q, constructed in the following way: given 
(n, k) G An 

{X n . k {u)) e ife^e/s+i, 

1 if Un,k+i <Pe, e = e k+1 and X nik (u) G K e , 

1 if u n ,k+i <Pe, e = e k+1 and X„ )fe (u) £ if e , 

otherwise. 

Furthermore notice that for all n there exists the limit X„ k f° r going to infinity. 
We construct X n+ i t o as this limit. We call such a process an Fif^-process. 

All DLR random cluster measures are invariant under this process. Indeed if 
4> G TZp.q and X n ,k has law (f> then X n ,k+i has the same law. Hence, for h > k and 
m = n, or for m > n, X m ^ has law </>. 

Finally, for any integer N < and any u 6 Jl we denote with ^" A i A '' 1 the 
F-fT^g-process constructed starting with Xn,o( u ) = w for u £ W. 

In the main result of this paper, under suitable assumption on the parameters 
p, q and on the graph G, given a finite F C E and u G U we will show how to 
determine a integer N such that (Ag^ ) ' A '' l (u)) e does not depend on uj G £1 for all 
e G F. Moreover we show how to determine a finite region F' D F such that 
(XQ^ ) ' Ar ' ) (u)) e does not depend on the values of u n , k for e k £ F' . 



(X nik+1 (u)) t 
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Since the 7?. Pig -measures are invariant in this way we prove that they are all equal 
and we give a perfect simulation of them on any finite subset of E. In particular 
we prove 1Z Piq = W p , q = 

4. SlMPLICIAL GRAPH 

In this section we define the notion of simplicial graph and we prove some geo- 
metric properties of these graphs. It is possible that these results are already known, 
maybe with different notations, however we could not find a reference. 

The prototypical graph we have in mind is the graph h d whose vertices's are the 
elements of Z d and whose edges are the segments of length one joining them. More 
in general a simplicial graph will be the graph obtained by considering vertices's 
and segment of a polyhedral tessellation of M. d . 

Before giving the details we explain roughly the problem we want to consider in 
the next sections. Let G = (V, E) be a graph and color each edge of G white or 
back in a random way. Now consider a finite subset F of E. We want to determine 
a region H containing F such that for all e G F if there is a path of black edges 
joining the two end vertexes of e then there is a path contained in H of black edges 
joining the two end vertexes of e (see Theorem 4). Moreover we want H to be small 
as possible so that if the probability of an edge to be black is high then if F is finite 
then H is also almost surely finite (see Theorem 7). We construct first a set H by 
adding inductively to F white edges until its "boundary" is entirely composed by 
black edges and then we set H to be equal to the union of H with its boundary. 
The idea is that, in this way, if two points are in the boundary of H , then they can 
be connected by a path contained in the boundary and in particular of black edges 
(for the precise statement see Proposition 3). In this way when a path of black 
edges has its end vertexes in F can be replaced by a path of black edges with the 
same end vertexes and contained in H by replacing the pieces outside H with paths 
along the boundary. A first try could be to construct H by adding edges to F until 
there are white edges connected by a vertex to the set. In this way H would be the 
union of F with the connected components of the subgraph of white edges having 
non trivial intersection with F. However it is immediate to see that this H has not 
the required properties. For a simplicial graph we can construct H replacing the 
condition "to have a vertex in common" with a different condition. We construct 
H by adding to F white edges until there are white edges in the "boundary" of this 
set. The correct notion of "boundary" is defined in Section 4.3, now we explain it 
in the case of the simplicial graph h d . We say that an edge e' is in the "boundary" 
of an edge e if they are different and if there exists a ci-dimensional hypercube of 
side 1 with vertexes in 1 d containing both e and e' . In this section we prove that 
this notion of boundary has the required geometrical properties (see Proposition 
3). 

4.1. Definitions. Let A be a closed convex subset of R e which is the intersection 
of a finite number of closed half-spaces. Such a set will be called a convex cell and 
will be the starting point of our constructions. For such a set we can identify the 
subset of vertexes, edges and z-dimensional faces and we denote by Ai the set of 
i-dimensional faces of A. 

In some constructions will be useful to have a more general notion of cell. If A 
is a convex cell of dimension m and tp : A — > M. d is a piecewise affine continuous 
injective map we call the image a of A a cell of R d of dimension m. We define also 
the collections Uj = {<p(B) : B e Ai} and the datum of <7o, . . . , <J m will be called a 
■polyhedron. 
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The assumption piecewise affine on ip could be highly relaxed, however this 
assumption makes some of the arguments below more elementary and does not 
change the generality of our applications. 

A polytope V in R d is a collection of cells in M. d which intersect properly. More 
precisely V is the datum of sets Po, Pi, ■ ■ ■ , P m such that 

i) the elements of Pi are i-dimensional cells of R d ; 

ii) Pi is locally finite: this means that for all bounded regions R of R d aHR = 

for all a e P- L but a finite number; 
Hi) if we denote by Pj(a) = {t £ Pj : r C <r}, for all a 6 Pj and all j < i, then 
the collection 

V{a) = {P {a),...,Pi{a)} 

is a polyhedron and the set Pi(cr) will be called the set of i-faces of er; 
iv) for all a e Pi and r £ P, the intersection a n r is either empty or a union 
of faces of er and r. 

The elements of Pi will be called the i-cells of V and in particular we will call 
Po (resp. Pi) the set of vertexes (resp. of edges) of V. The union of all the cells of 
V will be called the support of V and will be denoted by supp (P). 

The graph G(V) = (V, E) associated to P is defined as follows: V = P , E = Pi 
and the end vertexes of an edge e 6 Bis the pair of vertexes contained in e. Notice 
that if x and y are two vertexes of P then they are in the same connected component 
of supp (P) if and only if they are in the same connected component of G(P). 

The simplest possible convex cell are the simplexes. The I dimensional standard 
simplex is the set S = {(xq, . . . , xg) E M. i+1 : Xi > for all i and J2i x i = !}■ 
A ^-dimensional convex simplex (resp. an ^-dimensional simplex) is the image 
of S under an affine (resp. piecewise affine and continuous) injective map. We 
notice that every cell can be obtained as the support of a polytope whose cells are 
simplexes. 

A refinement of a polytope P is a polytope V such that supp (P) = supp (V) 
and each cell of V is contained in a cell of P. 

4.2. Internal and external part of codimension one smooth polytopes. We 

say that a polytope C in R d is smooth if its support is smooth as a topological variety. 
In this case this means that for every x <G supp (C) there exists a natural number j, a 
neighborhood W of x in R d an open ball W of K d and a piecewise affine continuous 
map tp : W — > W which defines an homeomorphism between W and W, such 
that ip(x) = and ip(W n supp{C)) = {{h, . . . , t n ) e W : h = ■ ■ ■ = tj = 0}. If, 
moreover all the connected components of supp (C) have dimension n — 1 we say 
that it is a smooth polytope of codimension one. Notice that if C is smooth of 
codimension one in R d then for all x in supp (C) there exists a neighborhood W of 
x in M. d which is divided by supp (C) into two open connected components. Now 
we give a more global construction of these components. 

Let C be a codimension one smooth polytope in R d with a finite number of 
vertexes and let U be the complement of supp(C). For all x £ U we consider the 
set S X (C) of half-lines i starting in x and whose intersection with supp (C) is generic. 
More precisely we require that for all cells CTofCif^n<7^0 then a 6 Pj-i and 
t H a is a finite set and moreover this intersection is contained in the set of points 
of a that are linearly smooth: for all y G i n a there exists a neighborhood U of y 
and an hyperplane such that a n U = U OH. Then the parity of the cardinality of 
I n supp (C) does not depend on i 6 S X (C). Moreover this parity is locally constant 
on x € U. 

Hence we define the internal part of C, that we will denote by IntC as the set 
of points x e U such that this cardinality is odd and the external part, that we 
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will denote by EstC as the set of point such that this cardinality is even. Notice 
that IntC and EstC are two open subsets of R d with boundary equal to supp (C). 
In particular for all path joining an element of Int C with an element of Est C the 
intersection of 7 with supp (C) is not empty. Finally notice that IntC is bounded. 

Notice that for all x £ supp (C) if W is a neighborhood of x as in the beginning of 
this section, then the two connected components of W\supp (C) are the intersection 
of W with IntC and EstC. We say that a path 7 cross C in x if 7 (to) = x f° r some 
t and there exists sequences {t n } and {s n } going to t such that j(s n ) £ EstC and 
j(t n ) £ IntC for all n. 

Lemma 1. Let C be a codimension one smooth polytope in R d with a finite number 
of vertexes. Then for all x, y £ supp (C) such that there exists a path in Int C 
joining x and y, and a path in Est C joining x and y, then x and y are in the same 
connected component of supp (C). In particular if x,y are vertexes then they are in 
the same connected component ofG(C). 

Proof. Let a (resp. ft) be a path joining x and y in IntC (resp. EstC). Consider the 
connected component D of supp (C) containing x and let V be the polytope whose 
cells are the cells of C contained in O. T> is a connected codimension one smooth 
polytope with a finite number of vertexes. Notice that, by what noticed above, in 
a small neighborhood W of D we have that W\D has two connected components 
Wi and W 2 and we have W 1 = W n IntC and W 2 = W n EstC. Similarly W x and 
W2 must be the intersection with Int V and Est V, and both the possibilities 



can occur. Since a and ft do never cross C they also never cross V and we have 
that a is contained in Int£> and ft is contained in Est£> or the opposite. 

Hence the final point y of a and ft belongs to Int T> n Est V hence it is in D. □ 

4.3. Polyhedral tessellation and simplicial graph. We say that a polytope 
V is a polyhedral tessellation of R d if supp (V) = M. d . In this case we say that 
G = G{V) is a simplicial graph. 

If A C Po we define the boundary A-p (A) of A as the set of vertexes v £ Pq \ A 
for which there exists, w £ A and a cell a of V, such that v,w £ a and we define 
G(A) as the graph whose set of vertexes is equal to A and whose edges are given 
by the element in Pi joining two elements of A. Notice also that if x £ A and 
y £ Po\ A any path in G(V) joining x and y intersects A-p (A). 

Proposition 2. Let V be a polyhedral tessellation of R d and let V be its set of 

vertexes. Let A be a finite subset of V and set B = V \ A. Let x,y £ A adjacent 
respectively to x', y' £ B. Assume now that x, y are connected in G(A) and that 
x' , y' are connected in G(B). Then x',y' are connected in G(A-p(A)). 

Proof. The proof follows exactly the same lines in the case of a tessellation using 
convex cells and in the case of a general tessellation. However in the first case all 
construction are more intuitive and direct. For this reason we give first the proof 
in the case of a tessellation made of convex cells and then we briefly explain how 
to change the proof in the general case. 

In the first step of the proof we assume also that all cells are convex simplexes. 
We construct a smooth polytope C of codimension one that separate A and B in 
the following way. 

We define the set of vertexes Co of C as the set of middle points of edges joining 
an element of A and an element of B. For all i dimensional simplexes a of C which 
contain an element of Co the convex envelope of er n Co is a i — 1 cell. Let Cj_i 
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be the collection of these cells and let C be the polytope whose i dimensional cells 
are given by C, for i = 0, . . . , n — 1. By construction C is a smooth polytope of 
codimension one with a finite number of vertexes. Notice also that a path in G(A) 
or in G(B) will never cross supp(C). 

Let now u,v G Co be the middle points of the edges joining x, x' and y and 
y' respectively. By Proposition 1 u and v are in the same connected component 
of G(C). Hence there exists a sequence of vertexes Wq — u, W\, . . . , w m = v in Co 
determining the path connecting u and v in G(C). Furthermore let ^ £ i and 
t- e 5 be such that Wi is the middle point of the edge joining ij and ^. Then 
t' = x' , t'i, . . . , t' m — y' determine a path in A-p(A) joining x' and y' . 

Let now V be any tessellation with convex cells. We construct a sequence 
of refinements of V and of finite subsets of the vertexes of pW. 

• is P and A« = A. 

• is the tessellation obtained by adding a vertex v a in the barycentre of 
all 2-dimensional faces a G P2 which are not a simplex and adding the edges 
joining v a with the vertexes of a. Finally we set A^ = Au{vcr : <jC\A ^ 0}. 

• more generally given P^ x \ pM will be the tessellation obtained adding 
a vertex v a in the barycentre of every i-dimensional cell a of which 
is not a simplex and adding all the j-dimensional cells obtained by joining 
this vertex with the j — 1 cells contained in a. Finally we set A^ = 
A^-V U {v a lanAf'- 1 ) ^0}. 

Set V = V {d) and A' = A^ , and B' is the complement of A 1 in the set of vertexes 
of V'. Notice that all cells of C are convex simplexes and that A-p(A) = A-pi(A'). 
Hence we can apply to this situation what we have already proved. 

In the case of a tessellation whose cells are not convex some of the constructions 
we have described have to be modified. Indeed it does not make sense to consider 
the middle point of an edge or the convex envelope of the middle points in the first 
part of the proof or the barycentre in the second part of the proof. As an example 
we show which changes are necessary in the construction of the sequence V^ l K 

Suppose we have already constructed "P^ -1 -* and that all j-dimensional cells 
of P^" 1 ) are simplexes for j < i — 1. Now consider a i-dimcnsional cell a of 
p(* _1 ) which is not a simplex. Let tp : A — > a be the piecewise affine map which 
parametrise a where A is a convex cells. Notice that by induction all faces of A 
are convex simplexes. Now let v be the barycentre of A and that we can divide A 
into simplexes A T joining v with the faces of A. We obtain the new tessellation by 
considering the restriction of (p to the simplexes A T . □ 

The result we will use in the next section is a variation of Proposition 2 where 
the set of vertexes A is replaced by a set of edges. 

If H C Pi we define the boundary A-p (H ) of H as the set of edges e G P\\H 
for which there exists, d G H and a G P n such that e,d C o. Let also Vh = {v G 
P : v G e for some e G H} and let G{H) be the subgraph (V H ,H) of G(P). 

Proposition 3. Let P be a polyhedral tessellation of R d and H C P\ be a finite 
set. Let x,y G Vh be connected in G{H) and in G{P\ \ H). Then x and y are 
connected in G(A-p(H)). 

Proof. Let e\, . . . , e m be a path in H which join x and y and let fi, ■ ■ ■ , f n be a 
path in Pi \ H which joins x and y. Let x = uq, . . . , u m = y and x = Vo, . . . , v n = y 
be the sequences of adjacent vertexes determined respectively by the path in H and 
in Pi \ H. Dividing the path into smaller pieces we can assume that Ui ^ vj for 
1 < i < m — 1 and 1 < j < n — 1 . 
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We construct a refinement Q of V replacing each vertex w with a small convex 
cell so that we will be able to apply Proposition 2 to this situation. 

More precisely for each vertex to we consider a "small" ball B w with centre w 
which is "far" from every cell of V not containing to. We set also D w to be the 
intersection of the boundary of B w with the edges in Pi containing w and A w to 
be the intersection of the boundary of B w with the edges in H containing w. We 
set 

<9o = ^oU (J D v and A = {u t : i = 0, . . . , m} U (J A w . 

wGP w£P 

We now describe the higher dimensional cells of Q. Let F w be convex envelope of 
the set D w and let M n be its boundary and N w its open part. If B w is small enough 
for every t-dimensional cell a of V containing w the intersections o~ w — a n F w and 
erf = a ("1 M ro are respectively an t'-dimensional cell and an (i — I)-dimensional 
cell. Finally for all z-dimensional cells a of V the set a = a \ \J W N w is also an 
^-dimensional cell. The tessellation Q is the collection of the cells a, a w and erf 
for a a cell of V and to G Po- 

Notice that if a is an edge in G(Aq(A)) then a is an edge in G(Ap(H)). 

We now apply Proposition 2. The sequence of edges 

(^1 )uo 5 ^1 ) (^1 )ui ) (^2)^11 ^-2i • • • 1 &m 7 {^m^)u m 

is a path in G(A) joining x and y. Let x' = and t/' = (f n )y then they are two 

vertexes adjacent respectively to x and y and the sequence of edges 

fl, (fl)v! j {f2) Vl , /2, • • • , (/n)%_n/n 

is a path in G(Qo \ -4) joining x' and j/'. Hence there exists a path 

Ti,...,T ri ,CTi,T ri+ i,...,T r2 ,Cr 2 ,...,T r<j 

in G(Aq(A)) joining x' and y' where we assume that r, is of the form r\ w or 77^ for 
all i. Then cri, . . . , ct^-i is a path in G(A-p(H)) joining x and y. □ 

5. Coupling of the random cluster measure at low or high 

temperatures 

In this section we fix a polyhedral tessellation V in M. d and we denote by G = 
(V,E) the underlying simplicial graph. We introduce a new graph G* = (V*,E*) 
where V* = Z x V and E*=ZxEUZxl/ where if x, y are the end vertexes of 
e e E then (n, x), (n, y) are the end vertexes of (n, e) G -E* and if v G V and n G Z 
then we denote with e n> „ the corresponding element in E* and its end vertexes 
are (n, u), (n + 1, w) G V*. Notice that G* is the simplicial graph of a polyhedral 
tessellation V* of = R x K'' whose z-dimcnsional cells are the collection of 

the cells of the form {n} x a where a is an z-dimensional cell of V and n G Z and 
of the cells of the form [n, n + 1] x r and r is an (z — l)-dimensional cell of V and 
n G Z. We consider G as the subgraph G x {—1} of G*. 

Fix u eU. We define the following coloring of the edges E* of G*: 

W/ = W(u) = {(n, e fe ) G E 1 * : (n, fc) G ^4 and p efc < u„ :fe < 1}, 

M = M{u) = {(n,e k ) G E* : (n,k) G A and p efc < u„ :fe < p e J, 

B = =E*\ (M(u) U VF(u)). 

We define also B n = {e G E : (n, e) G B}. We say that the elements of B are black, 
the elements of M are gray and the elements of W are white. 

Given a subset F of M(u) U W(u) we define the cluster of white or gray edges 
Cp = Cp (u) as the minimum subset H of E* containing F such that if e G E* and 
e G A-p* (if ) then e G More in general if F C £7 is not necessarily a subset 

of M U W we define G$ = F U G™, B . Equivalently if F C M we can construct Cf 
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inductively by adding the white or gray edges "near" to F as follows. Let D = F 
and A+i = A U {A v , (Di) \ B) then Cf = \J l D t . 

Notice that as the numbers p e grow the probability that Cp is finite increase. As 
we will see in the last section in the Ising model the p e 's are related to a parameter 
called temperature and when this parameter is small the p e are closer to 1. For this 
reason we refer to the case in which Cp is finite as the situation at low temperature. 

Assume now that Cp is finite. In this case we define N w (u,F) as the biggest 
negative integer N such that Cp cEx [N +1,-1}. We define 

H% = H™{u, F) = {eeE: (n, e) e C^ or (n - 1, e) e Cp 1 } (5) 

and we set H%(u, F) = H™ U A V (H%) and finally H w (u, F) = \J n {n} x H%. 

Theorem 4. Fix u eU such that Cp(u) is finite and an integer N < N w (u,F). 
Let v! G Un such that u' n k = u n ^ k for all (n, e k ) G H w (u, F). Then for all ui, ui' G f2 
we have 

for all e G H^(u,F). 

Proof. Let r\ nM = X^f\u) and n' nk = X^' k ' N) (u'). Let also C = Cf and C n = 
{e G E : (n, e) G C}, H n = C n U C„_i and H n = H n U A P (H n ). We notice first 
that A-p(H n ) C B„n B n _\. Indeed let e G A-p(H n ), then there exists e' G i?„ 
and a cell cr in V such that e and e 1 are contained in a. If e ^ _B„ then by the 
definition of C and of the cells of V* we have that (n, e) G C, hence e G C„ which 
is in contradiction with e G A-p(H n ). Similarly we get an absurd if e ^ B n _\. By 
the definition of a Fif^-process this implies that 

Ap^jc^jn^) (6) 

for all to > N and ft, > 0. 
We will prove that 

(Vm,h) e - (»7^,fc) e ( 7 ) 

for all (m, ft) G ^Iat and for all e G ff m . We prove this by induction starting with 
to = N and ft = 0. For m < N w (u,F) and for all ft the equality (7) is trivially 
satisfied since C m = 0. 

Now we prove that if equation (7) holds for m = N, . . . , n — 1 and for all h and 
for to = n and ft = 0, 1, . . . , k then it holds also m — n and h = k + 1. Let e G 
If e ^ ek+i then by induction and definition of FKp q -process we get 

{Vn,k+l) e = {Vn,k) e = Wn,k) e = {v'n,k+l) e 

proving the claim. If e = e^+i and e G H n we compute {rj n] k+i) ek i - If u n ,k+i > 
Pe k+1 then wc have (rj n k+i) — and similarly for n' . If u n k+i < Pe k+1 then we 
have {Vn,k+i) e = 1 and similarly for r/' . If p efc+1 < u n . k+1 < p ek+1 we need to 
prove that 

Vn,k&K ek+1 iff n' n k G ife H i • (8) 
Let x and y be the end vertexes of e k +i and assume that there is a path 7 : E\, . . . , e m 
in E(-q n ^ k ) \ {e k+1 } joining x and y. 

If the path is contained in H n — C n U C„_i U A-p(H n ) then wc prove that the 
same path is contained in E{n' n k ) \ {e k+ i}. Let e r be an edge of the path. By 

induction {r] n .k)e r = {v'n,k)e r hence e r G E (v' n ,k)- 

If 7 is not contained in _ff„ wc show there is another path joining x and y 
contained in H n . By (6) and the fact that p ek+1 < Un,k+i < Pe k+1 we have e k +\ G 
_ff„. Let D be the connected component of H n containing e k +i- Let s\, . . . , £,_i G 
D, Ei, . . . ,sj £ D and £j + i G £>. Let x' be the vertex common to £,_i and £j and y' 



PERFECT SIMULATION FOR THE RANDOM CLUSTER MODEL 



11 



the vertex common to ej and £j+i- We can apply Proposition 3 and we construct a 
path /3 in A V (D) C A v (H n ) joining x' and y'. Since A v (H n ) C E(r) n , k ) n F« i/s ) 
we can replace 7 with the path 7' : £1, . . . , £,_i, fi, £j+i, ■ ■ ■ , £ m - Repeating this 
process we see that we can substitute the path 7 with a path entirely contained in 
D U A-p (D) C H n as claimed. Hence we are reduced to the previous case. 

Finally we prove that if (7) holds for a fixed m and all h > 1 then it holds also for 
m + 1,0. Let e G C m +i- If e G H m this follows by definition of the FAT^ -process. 
If e = e r G Hm+i \ C m then e G F m otherwise e would be an element of C m ). Then 

0?m+l,o)e r = ( J 7m,r)e r = 1 = (v'm,r)e r = (?7m+l,o)e r 

proving the claim. □ 



5.1. Coupling at high temperatures. We give now a similar result correspond- 
ing, in the Ising model, to high temperatures. 

Let G — (V, E) be a countable graph (in this case we do not assume simplicial). 
Define G as the graph with set of vertexes V — Z <0 x V and edges V — Z <0 x E 
where if e G E has end vertexes x, y then the edge (n, e) has end vertexes (n, x) 
and (n,y). We consider G as the subgraph G x { — 1} of G* . For all u G W define 
M(u), M^(u) as in the previous section. 

For a subset H of i? we denote by the set of edges which are not in H and 

which have a vertex in common with an edge in H. 

Fix u eM and a subset F E* . If F C M(u) we define i/ie cluster of black or 
gray edges C h F = C b F (u) as the smallest set C of E containing F and such that for 
all (n, e) G E\C if either (n — 1, e) or (n, e) or (n + 1, e) have a vertex in common 
with an edge in C then (n, e) G W{u). If F is not necessarily contained in M(u) 
we define C b F {u) = G FnM ^(u) U F. For all n < we define 

H b n = H b n {u,F) = {e G E : (n, e) G C£(u) or (n- l,e) G C7^(w)} 

and we set H^(u,F) = H b UT(H b ) and H b (u,F) = {J n {n} x H*. 

Finally if C b F {u) is finite we define N b (u,F) as the biggest negative integer N 
such that C£ C E x [TV + 1, -1]. 

Theorem 5. Fix u Gli such that C F (u) is finite and an integer N < N b (u,F). If 
e and u' G Un is such that u' n k — u n ^k for all (e/s,n) G H b (u, F) then 

for all e G Ht^F). 

Proof. The proof follows exactly the same strategy of the proof of Theorem 4. 
However in this case it is simpler since we do not have to use the result of Sec- 
tion 4. We give here only the main lines. Indeed an argument analogous to proof 
of equation (6) gives 

T(H b (u, F)) G W{u). (9) 

Then we prove the equality (7) by induction as in the proof of Theorem 4. Also 
in this case we are reduced easily to prove the equivalence (8). This equivalence 
is easier in this case, since, by (9), we have that E(r] nt h) C H b for all n > N and 
similarly for r/' so we can assume that the path joining the extremal point of e k +i is 
contained in Fl\ without using any further result, while in the proof of Theorem 4 
we need Proposition 3. The remaining argument are completely similar to the proof 
of Theorem 4. □ 
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6. Assumptions for the finiteness of clusters 

In this section, given a graph G, we present some conditions on it and on the 
parameters p such that the cluster Cp of Theorem 4 is almost surely finite or such 
that the cluster C F of Theorem 5 is almost surely finite. 

We start by recalling a general Lemma. Let G = (V, E) be a graph and let 
7T = (n v ) ve y be an element of [0, l] v . Consider the product measure on the space 
f^G = {0, 1} V such that P(u> v = 0) = n v . For each oj G fie let G[w] be the subgraph 
of G with set of vertexes V[w] = {v G V : w v =0} and with edges the set E[oj] 
of the elements of E joining two elements in V[w]. Moreover for each uj G SIq and 
for any i> G V set C^w] the connected component of G[w] containing v (possibly 
empty if v $ Y[u>}) and if n is a natural number let V„ in [w] be the set of vertexes 
in G v [co] whose graph-distance from v in the graph G v [u] is equal to n. 

Lemma 6. Let G and it be as above. For each vertex v G V let A v be the set of 
vertexes adjacent to v and set g v — J2veA n v- If g — snp{g v : v G V} < 1 then 
P({uj G n G : V„, n [w] ^ 0}) < g n . 

Proof. Define the random variables Z Vin [u] = card (V^nfu;]) . The conditional mean 
value E(Z v ^ n+ i\Z Vin ) verifies E(Z v ^ n+1 \Z v ^ n ) < g Z ViU , in particular the sequence 
{Z Vtn } n is a supermartingale. Hence the mean value E(Z v ^ n ) is less or equal to g n . 
By Markov inequality the claim follows. □ 

6.1. Assumptions for the finiteness of clusters at low temperatures. Now 

we give conditions on G and p such that the cluster Cp (u) defined in Section 5 is 
finite for almost all u eW. We fix a polyhedral tessellation V of R d . Let G = (V, E) 
be the associated simplicial graph and let V* and G* = (V*,E*) be defined as in 
Section 5. 

For the proof of our next Theorem we introduce a new graph G = (V, E) defined 
as follows: V — E* and two elements e,e' G V are joined by an edge in E if and 
only if there exist a cell in V* which contains e and e'. 

If H C E* and u e U we define also the subgraph G(H, u) = (V(H, u),E(H, u)) 
of G whose set of vertexes are equal to V(H, u) = H\B{u) and whose edges E(iJ, u) 
are all the edges of E joining two vertexes in Y(H, u). 

For e e E define 

g e = 2(i-p e ) + 3 Y, (i-iv), (io) 

e' G A-p ({e}) 

where p e is defined in (2). 

Theorem 7. Let G = (V,E) be a simplicial graph in R d . Assume that p e > for 
all e e E and that 

Um sup § e < 1. 

A T-fc e^A 

Then for all e G E the cluster C™{u) is finite for almost all u^U. 

Proof. First we do a preliminary remark: the event / = \J e ^p{u G IA : card (C™(w)) = 
oo}, is in the tail cr-algebra. Therefore, by Kolmogorov 0-1 law, this event has 
probability zero or one, in particular to prove our claim it is enough to prove that 
P(I) < 1. 

Define F n — {e G E : g e > 1 — i}. By assumption there exists n such that 
card(i^ no ) < oo and set F = F 7lQ and g = 1 - Set also F = F U Aj>(F ) and, for 
£>1, define F = ZxFcZxEcE*. Fix an edge eeE and, for £ > 1 set 

= [-1,-1] x B G (e,£) cZx£ 

and set also Se — card(A£). 
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Choose £o such that 3 cardjF) ^g < 1 an( ^ fo r £ > £ + 1 ; define the events 

Wg = {ueU : u„,fc < p efc for any (n, e fc ) £ [-f - - 5 £ + £ ] x F} 

Xe = {u €U : card (C£ («)) = oo}. 

We notice that the probability pj = P (VV«) docs not depend on £ and it is a 
positive constant being p e > for any e <G £7. We notice also that C and 
that their union is equal to I. Therefore 

P{I) = P^^=VmP(X l ). 

We also define the events 

ye = {u G U : there exists a sequence ei, . . . , e m G W(u) U M(u) \ F 
such that ei is adjacent to e i+1 in the graph G, 
ei G A £ and e m ^ A £+S J 
Ze.i = {u eW : there exists a sequence ei, . . . , e m G W(u) U M(u) such that 
ei is adjacent to e i+i in the graph G, 

ei G {-£ - S*^ + z} x F, e m ^ A^ +Sf and e» ^ F for 1 < i < m}. 
Finally define Zg t i = Zg t i and Zgj = Zgj \ -2^,-1 for i > 1. It is clear that 

c % U (J Z e ,i, 

i=l 

in particular P(^|W,) < P(^|W,) + Eti ^C^W/). 

Now we notice that the events Wg is decreasing meaning that if u G Xg and 
u' G W is such that h < u n ,h for all n, /i then u' G Wg. With a similar definition 
the events yg and Zg^ are increasing. Hence, by the FKG inequality we obtain 
P{ye\We) < Ptye) and P(Zg ti \Wg) < P(Z e>i ) (see [Gri99], Chapter 2). Hence, 
noticing that P(Zg. t \Wg) = for i < £ , we get P(X e \W e ) < PQ>/) + £ j=/o 

Now we estimate P(3^) and P(Zg^) using Lemma 6. We start with yg. Consider 
the random graph G(F c ,w) and set 7r e = 1 — p e . Notice that yg G {J eeAf {u G U : 
Y(F c ,u) e ,s e ^ 0} hence using Lemma 6 we get 

P( k yt)<ca I d{A e )g St = S e g St , 

for ^ large enough. For Zg.i we proceed in a similar way. Consider again the 
random graph G(F c ,w). Notice that in the sequence e\,...,e m which appears 
in the definition of Zg^ the subsequence e2-..,e m _i is in F c and e2 G Fg.i := 
{-£ - Sg + i - 1, -I - S e + i, -£ - Sg + i + 1} x F. Hence Z e<i C U eG F t .{u <E U : 
V(P C , w) ej i ^ 0} and using Lemma 6 we get 



P{Zg,i) < card(F^)^ = 3 card(P) 5 l 
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Recall that Pq — P(We) does not depend on £, hence we have 

p(x t ) = p(Xi\wt)p{m) + p{Xi\w c e )p{w c t ) < p{wi) + p(Xt\m)P(Wt) 

oo 

<i-p + p (p(y e ) + J2 p ( 2 i,i)) 

i=io 

<i-Po + Ms e9 ^ + 3 -^l g ^) 

l-g 

<l-^ + Se 9 s 'Po. 
Finally notice that lim^oo Sg g Se = 0. Hence 

p 

lim P(X e ) < 1 < 1 

as claimed. □ 

6.2. Assumptions for the finiteness of clusters at high temperature. Let 

G = (V, E) be a countable graph of finite degree. We remark that in this case we 
do not need to assume that G is a simplicial graph. For all e £ E define 

g e = 2p e + ?> Pe'- 

e'er({e}) 

Theorem 8. Let G be a countable graph of finite degree. If p e < 1 for all e 6 E 
and 

lim sup ^ < 1 

then for all e E E the set C^(u) is finite for almost all u eU. 

The proof follows exactly the same lines of the proof of Theorem 7, however we 
do not need any result from Section 4. 

7. Perfect simulation of the random cluster measure at low or high 

temperature 

As an application of the previous results we now explain how to prove uniqueness 
of the random cluster measure and how to obtain a perfect simulation of the random 
cluster measure using the results of the previous sections. We explain these results 
in the case of low temperatures. The case of high temperatures can be obtained in 
a similar way. In this section, from now on we assume that G is a simplicial graph, 
that p e > for all e E E and that limAfs sup e ^ A g e < 1. The uniqueness proved in 
the following Corollary is well known at least in the case of L d . 

Corollary 9. Assuming the hypotheses above the random cluster measure on G is 
unique. 

Proof. Let (f>, <j)' be two DLR random cluster measures. 

To prove that <j) an d 4>' are equal we prove that for each finite subset F of E 
the projections (f>F and 4>' F of <p and <j>' onto {0, 1} F are equal. We denote also by 
X { p' u) e {0, 1} F the projection of . 

Let co be a random variable with law (f> and u/ a random variable with law <p' . 
By Theorem 4 we have that 

\\xf^ ] - X^^Wtv <P(ueU: N w (u,F) < N) 

where the lefthandterm is the total variation distance between the law of X^'^ and 
X( F N '"'\ Recall now that as noticed in Section 3 the DLR random cluster measures 
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are invariant under a -FA"p ^-process. Hence, since to has law </>, the random variable 

Xp?'^ has law <j>p and Xp 1 '" -* has law <jjf F . Hence we get \\<j>F — 4>' f \\tv < P(u e 
U : N w (u,F) < N). Finally by Theorem 7, Cf{u) is finite for almost all u G U. 
Hence P(u G U : N w (u, F) < N) goes to zero as N goes to infinity. Hence <j> and 
4>' are equal. □ 

Now, under the same assumptions, given a finite subset F of E we briefly describe 
an algorithm which furnishes a sampling of the random cluster measure on F. 

Let be given a generator of independent random numbers u n ,fc. The algorithm 
takes F as an input, a suitable description of the tessellation V, and gives as output 
subsets C, H of E* , a subset of E and a configuration Y G {0, l}^- 1 . It uses 
also local variables D,D',F',L and generates u n ^ for (n, e^) G if. We uses also 
the notation of B(u), M(u),W(u) for black, gray and white edges introduced at the 
beginning of Section 5. We describe the algorithm with the following pseudocode. 

Step 1: generate the random numbers u n ^ for n = — 1 and e k G F and set 

F' = {-1} x F; 
Step 2: set D = F'\B(u)); 

Step 3: set L = DliA-p, (D) and generates the random numbers u n u for (n, e^) G 
L\(DUF')- 

Step 4: set D' = D U (L \ B(u)) as in the construction of the sets Di at the 

beginning of Section 5; 
Step 5: if D' ^ D assigning to D the value given by D' and goes to step 3; 
Step 6: if D' = D then C = D UF' and compute N = N w (u, F) as in Section 5; 
Step 7: use the formula (5) to compute the sets H™ and define H = {J~= N {n} x 

H n where H n = H% U A V (H%) as in Section 5; 
Step 8: generate the random numbers u n ^ for (n, ek) & H \ C; 
Step 9: e G H-i use the process described in Section 3 to compute the value 
Y e = X^^^ where uj e = for all e G Hn. As explained in the proof of 
Theorem 4, by Proposition 3, for this computation it is enough to know 
the value of u only inside the region H. 
Notice that under our assumption, almost surely, this algorithm will end in a 
finite number of steps. Notice also that the output C is the set C F (u,F), H is 
the set H w (u, F). Finally as explained in proof of uniqueness above Y has law <f>F 
where 4>f is the projection onto {0, 1} F of the unique random cluster measure on 
{0,1}* 

In the case sup eg£; g e < 1 it can be easily proved that the average complexity of 
this algorithm goes linearly with the cardinality of F. 

8. Applications to the Ising and Potts Model at low or high 

temperatures 

We apply the results obtained in the previous sections for the random cluster 
measure to the construction of a perfect simulation of the Ising and of the Potts 
model with free boundary condition. 

Let G = (V, E) be a countable graph of finite degree. We briefly recall how the 
Ising model with free boundary conditions can be obtained from the random cluster 
measure with q = 2, see the original paper by Fortuin and Kasteleyn [FK72], the 
book of Grimmett [Gri04], Chapter 1 or the book of Newman [New97], Chapter 
3 for the details. Let (3 G (0, oo) be the parameter of the Ising model called the 
inverse of the temperature and let J e for e G E be the positive parameters defining 
the Ising model called the interactions. In the random cluster measure choose q = 2 
and p e = 1 — e^ 13,7 " for all e G E. Now we color, with + or — , each vertex of G 
using the following rule. Consider the subgraph G{lo) = (V, E{uj)) of G where to is 
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selected using the random cluster measure 4>^ 2 an d E(w) is defined as in Section 
2. Color the connected components of G(w) with + or — with an independent 
and uniform probability 1/2. Finally color each vertex with the same color of the 
connected component containing the vertex. A similar construction can be used to 
obtain the Potts model with free boundary conditions with the only difference that 
there are n colors, and we have to use the random cluster measure (p® n . Obviously 
the Ising model is a particular case of the Potts model so that we will study this 
second one in what follows. 

We explain how to obtain a perfect simulation of the Potts model with free 
boundary condition in the case of low temperature. The case of high temperature 
is completely analogous. So assume that the hypotheses of Theorems 4 and 7 are 
satisfied: G is a simplicial graph associated to the tessellation V, p e > for all 
e G E and Wm^E sup e ^ g e < 1. Given a finite subset W of V we briefly describe 
an algorithm which furnishes a sampling of the Potts model on W. 

Let F' be the set of edges having at least one end vertex in W and let F be a 
connected set of edges containing F' . Now we use the algorithm explained in the 
previous section to obtain a sampling of the random cluster measure. In particular 
given u we produce a set H = and a configuration Y = Xqq N ^ (see 

Theorem 4 for the explanation of this notation). Define also K as the connected 
component of H containing F and K as K U A-p(K). Let Vg be the set of end 
vertexes of the edges in K. Let also GK(Y) be the subgraph of G(Y) having as 
vertexes the set Vg- and as edges the set K n E(Y). Notice that, as in the proof of 
Theorem 4, if e G A-p(H) then Y e = 1. Hence, by Proposition 3, if x,y G Vg- then 
they are in the same connected component of the graph G(Y) if and only if they are 
in same connected component of the graph GK(Y) which under our assumption, by 
Theorem 7 is almost surely a finite graph. We determine the connected components 
of GK{Y) and we color each connected component as prescribed by the Potts model. 
Finally we color each vertex in W with the same color of the connected component 
containing it. 

It is easy to translate this description of this algorithm in an actual pseudocode 
as we have done in the previous section. 
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